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Abstract 

We study r — > [JtViUi, i = e,fi,T decay in the model III version of the two Higgs 
doublet model. We calculated the BR at the order of the magnitude of 10~ 6 — 10 
for the intermediate values of the Yukawa couplings. Furthermore, we predict the upper 
limit of the coupling for the r — h°(A°) — t transition as ~ 0.3 in the case that the BR is 
~ 10~ 6 . We observe that the experimental result of the process under consideration can 
give comprehensive information about the physics beyond the standard model and the 
free parameters existing. 
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1 Introduction 



Lepton flavor violating (LFV) interactions are interesting since they do not exist in the standard 
model (SM) and give strong signal about the new physics beyond. Such decays have reached 
great interest at present and the experimental search has been improved. Among LFV decays 
the ones existing in the leptonic sector are clean theoretically in the sense that they are free 
from the nonperturbative effects. The processes /i — > ej, r — > e(/i)7, r — > eee, r — > e/i/i, are 
the examples of LFV interactions. There are on-going and planned experiments for \x — > 
(t — > fi'-f) and the current limits for their branching ratios {BR) are 1.2 x 1CT 11 [lj (1.1 x 1CT 6 
|2j). The numerical estimates predict that the BR of the proceses r — > eee, r — > ejifi are at the 
order of the magnitude of 10~ 6 [3J, which is a measurable value in the LEP experiments and r 
factories. 

In such decays, the assumption of the non-existence of Cabibbo-Kobayashi-Maskawa (CKM) 
type matrix in the leptonic sector forbids the charged Flavor Changing (FC) interactions and 
therefore, the physics beyond the SM plays the main role, where the general two Higgs doublet 
model (2HDM), so called model III, is one of the candidate. This model predicts that the 
LFV interactions can exist at loop level and the internal neutral higgs bosons ho and Aq play 
the main role. There are number of Yukawa couplings which describe the strength of the 
interactions lepton-lepton-neutral Higgs particle, appearing in the loops and their strength can 
be determined by the experimental data. In the literature, there are several studies on LFV 
interactions in different models. Such interactions are studied in a model independent way in 

in the framework of model III 2HDM [5,, in supersymmetric models |Ul 171 1%1 1§1 ITUl ITT1 112j . 

Our work is devoted to the study of the r — > fiViVi, i = e,fi,r decay in the model III 
version of 2HDM. This process exists with the help of internal scalar bosons h° and A to 
obtain the flavor changing transition r — > /i and the internal Z boson to get the output vv 
(see Fig. |TJ). The BR(r — > fiUiUi), (i = e, /i,r) is predicted at the order of the magnitude of 
10~ 6 — 10 -4 for the intermediate values of the Yukawa couplings, which are the free parameters 
of the model used and is strongly sensitive to the couplings for r — h°(A°) — r and r — h°(A°) —\i 
transitions. This can play an important role in the determination of the upper limits of these 
couplings, especially the one for the r — h°(A°) — r transition. Notice that there have been 
some experimental studies on this process in the literature |T5] . 

The paper is organized as follows: In Section 2, we present the theoretical expression for 
the decay width of the LFV decay r — >• fiViVi, i = e, fi, r, in the framework of the model III. 
Section 3 is devoted to discussion and our conclusions. 
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2 t — > fi V{ i>i decay in the general two Higgs doublet 
model 



Type III 2HDM permits flavor changing neutral currents (FCNC) at tree level. The Yukawa 
interaction for the leptonic sector in the model III is 

£y = r)FjkL<t>iEjR + tfjliLfaEjR + h.c. , (1) 

where i,j are family indices of leptons, L and R denote chiral projections L(R) = 1/2(1=1=75), 0j 
for i = 1,2, are the two scalar doublets, Ul and Ejr are lepton doublets and singlets respectively. 
Here <\>\ and 02 are chosen as 
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By considering the gauge and CP invariant Higgs potential which spontaneously breaks SU(2) x 
U(l) down to £7(1) as: 

V((f>i,(h) = c 1 (<J ) tfa-v 2 /2) 2 + c 2 (<J)tfo) 2 

+ +C 3 [(0+0! - V 2 /2) + <P+H 2 + C4[(0^l)(0 2 f 2 ) - (0^0 2 )(0 2 f 1 )] 

+ c 5 [Re(<Pt<p2)] 2 + c 6 [Im(<Pt<p2)] 2 + c 7 , (4) 

with constants q, i = 1, ...,7, H\ and H 2 are obtained as the mass eigenstates h° and A 
respectively, since no mixing occurs between two CP-even neutral bosons H° and h° in the 
tree level. Therefore it is possible to collect the SM particles in the first doublet and the new 
particles in the second one. The part which produce FCNC at tree level is 



£y,fc = ifjkL&EjR + h.c. . 



(5) 



Here the Yukawa matrices £^ have complex entries in general. Notice that,in the following, we 
replace ^ E with £j| where "N" denotes the word "neutral". 

Now, we consider the lepton flavor changing process r —>■ \i v v and we expect that the main 
contribution to this decay comes from the neutral Higgs bosons, namely, h and i n the loop 
level, in the leptonic sector of the model III, (see Fig. Q). The general effective vertex for the 
interaction of off-shell Z-boson with a fermionic current is obtained as 



T (REN) (r ^^ Z * ) = fl ^ + f 2 7m75 + h GiiX + h a ^ 



(6) 



where k is the momentum transfer, k 2 = (p — p') 2 , p (p') is the four momentum vector of 
incoming (outgoing) lepton. Taking into account all the masses of internal (m,i) and external 
leptons (mj 1 , mi 2 ), the explicit expressions for the functions /i, / 2 , / 3 and / 4 are 
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L l{2}, (A ) = (AO) + (™i - m h {2} (! - &)) 

L P r (A ) = m ft°(AO) (1 -a; - (x + y) - k 2 xy , 
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The parameters cy and are c A = — \ and cy = 4 — sin 2 6*vi/- In eq. © the flavor changing 
couplings represent the effective interaction between the internal lepton i, {i = e,fj,,r) 
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and outgoing (incoming) j = 1 (j = 2) one. Here we take only the r lepton in the internal line 
and we neglect all the Yukawa couplings except £^ tt and Cn,t/i m the loop contributions (see 
Discussion section). Notice that the parameter [i in eq. (JJJ) is the renormalization scale, the 
functions fi, / 2 , ^3, f± are finite and independent of \x. 

The matrix element M for the process r — > /i z/j z/j, i = e, //, r is calculated in the framework 
of the model III, by connecting the r — > fi transition and the z/^z/j output with the help of the 
internal Z boson. Finally the decay width V is obtained in the r lepton rest frame using the 
well known expression 

dT = | M | 2 s\ P - y Pl ) fr . ^ , (10) 

2m T 1 1 ^ ,y (2tt) 3 2^' 1 ; 

where p (pi, i=l,2,3) is four momentum vector of r lepton (/x lepton, incoming u, outgoing v). 

3 Discussion 

In the case of vanishing charged interactions, under the assumption that CKM type matrix in 
the leptonic sector does not exist, LFV interactions arise with the help the neutral Higgs bosons 
h° and A , in the framework of model III. In this scenario the Yukawa couplings = 
e, fi, t play the main role in the determination of the BR of the processes under consideration. 
Since these couplings are free parameters of the theory, they should be restricted by using the 
present experimental limits of physical quantities, such as BR of various leptonic decays and 
electric dipole moments (EDM), anomalous magnetic moments (AMM) of leptons. In general, 
these Yukawa couplings are complex and they ensure non-zero lepton EDM. 

Now, we briefly discuss the case that the known light neutrinos Vi as massive particles and 
the lepton numbers Li denote the leptons of i th family are not conserved. If this is so, the lepton 
sector is an exact analogy to the quark sector and there exists a similar CKM type matrix, 
Maki-Nakagawa-Sakata (MNS) matrix V[ v ^3], that its elements are measured in neutrino 
oscillation experiments. It has been shown that the mixing between the muon neutrino and 
the heaviest mass eigenstate of the neutrino sector, the V^3 element, is nearly maximal [To1ll6j . 
The experiments on solar neutrinos J To] . 17\ (the reactor experiments such as CHOOZ [TH] ) 
predicted the mixing between electron neutrino and the second heaviest mass eigenstate of the 
neutrino sector, the V e 2 element (the heaviest mass eigenstate of the neutrino sector, the V e ^ 
element). Notice that the corner element V e s is much smaller than the others. On the other 
hand some off-diagonal elements of MNS matrix, such as V^3 are large and the MNS matrix is 
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far from diagonal in contrast to the CKM matrix (see for example jT2] for more discussion on 
lepton mixing). 

With the inclusion of MNS matrix, the existence of the lepton flavor violating (LFV) pro- 
cesses, t — ► fi transition in the present work, in the SM, would be possible. However, we 
expect that the tiny masses of the internal neutrinos bring small contribution even with the 
choice of maximal mixing in the leptonic sector. (See for example [20]; for the discussion of the 
existence of the MNS matrix and its effects on a special LFV process). In any case, the process 
r — > fi Vi Ui, i = e, /i, r should be examined theoretically if the lepton mixing is switched on. 

In the framework of the model III the r — > \xZ* transition can be switched on with the 
internal neutral Higgs bosons h° and A , and internal leptons e,fi, r. This brings unknown 
free parameters ^N,u,j an d £,n,tji hJ = e ->J JL i T - By assuming that only the internal r lepton 
contribution is considerable, the Yukawa couplings which does not contain r index can be 
neglected. Such a choice respects the statement that the strength of the couplings are related 
with the masses of leptons denoted by the indices of them, similar to the Cheng-Sher scenario 
|2*T] . Furthermore, we take symmetric with respect to the indices i and j. Finally, we 
are left with the couplings £jv TT and £jf r „, which are complex in general. Notice that in the 
following we use the parametrization 

and we present the numerical values of some parameters used in the calculations as a table: 



m T 


1.78 (GeV) 


m z 


91 (GeV) 


m w 


80 (GeV) 




V0.23 


G F 


1.16637 x 10- 5 (GeV~ 2 ) 


T z 


2.49 (GeV) 


r r 


2.27 x 1(T 12 (GeV) 



Table 1: The numerical values of the physical parameters used in the numerical calculations. 



The measurement of the BRs of r — > \i Vi % = e, /i, r decays [IS] are based on counting the 
number of candidate jets and correcting for efficiency and event selection. In additions to this, 
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the backgrounds coming from taus decaying to hadrons or cosmic rays should be detected. In 
the process we study the output contains z/j z/ i; i = e, //, r and the extraction of this output from 
the most probable one z/ M v T (BR ~ 17.37 ± 0.09% |22|), which exist in the SM theoretically, is 
difficult from the experimental point of view. 

In this work we studied the BR of the process r — > fiViVi, i = e,fi,r and we used the 
numerical value ^/v-™ m the interval 5 (GeV) < |Cjv r J < 25 (GeV). Here, the upper limit for 
the coupling |£jv )T/ J has been estimated in [23| as ~ 30 GeV. In this work it is assumed that 
the new physics effects are of the order of the experimental uncertainty of muon AMM, namely 
10~ 9 by respecting , the new experimental world average announced at BNL [21] 

a M = 11659 203(8) x 10~ 10 , (12) 

which has about half of the uncertainty of previous measurements. Here we have not used any 
restriction for the coupling £jf TT except that we choose its numerical value larger compared to 
£jv r/1 . In addition to this, we expect the upper limit of £§ jTT by taking the upper limit of £^ T/i 
and the BR of the process as 10 -4 -10~ 6 . 

In Fig. El we present ^§ TT dependence of the BR for real couplings. Here solid (dashed, 
small dashed, dotted, dash-dotted) line represents the case for £§ tT „ = 5 GeV (10, 15, 20, 25 
GeV). This figure shows that the BR enhances with the increasing values of both couplings 
and it reaches to values at order of magnitude of 10~ 4 . Fig. |3] shows the possible values of 
TT and the ratio for the fixed values of the BR, BR = 10~ 4 (solid line) and BR = 10~ 6 

' ?JV,rr 

(dashed line). For = 0.1, the BR = 10" 4 is obtained if the coupling £§ TT ~ 150 GeV and 
the BR = 10~ 6 is obtained if the coupling £,§ TT ~ 50 GeV. The possible experimental search 
of the process r — > \i Vi i = e, r would ensure a strong clue in the prediction of the upper 
limit of the coupling £^ TT . 

Fig. H| represents the h° mass m h o dependence of the BR for the fixed values of Cn,th an d 
£§ >TT , ^§,t^ = lOGeV, £§ >TT = 100 GeV. This figure shows that the BR is sensitive to m^o and 
it decreases with the increasing values of m^a. 

Now, we take the coupling £§ TT complex 

£n,tt = \£n,tt\ e%dTT ' (13) 

and present the sin 9 TT dependence of the BR for £§ TT = 100 GeV for four different values of 
£jv !T „, namely £jf T „ = 5, 10, 15, 20 GeV (solid, dashed, small dashed, dotted lines) in Fig. El 
From this figure it can be shown that the BR is not sensitive to the complexity of the coupling 

FE 
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At this stage we would like to summarize our results: 

• We predict the BR at the order of the magnitude of 10~ 6 — 10~ 4 for the range of the 
couplings, £jv )TT ~ 30 — 100 GeV and £j| ~ 10 — 25 GeV. We predict the upper limit of 
the coupling for the r — h°(A°) — r transition as ~ 0.3 in the case that the BR is ~ 10 -6 . 
With the possible experimental measurement of the process t —> [iViUi, i = e, /jl, t, a 
strong clue in the prediction of the upper limit of the coupling £jf rr would be obtained. 
Notice that the upper limit of the coupling £^ rAl could be predicted previously (see for 
example [23)- This analysis also ensures a hint for the physics beyond the SM. 

• We observe that the BR is sensitive to the neutral Higgs masses m^o and m^o 

• We observe that the BR is not sensitive to the possible complexity of the Yukawa cou- 
plings. 

Therefore, the future theoretical and experimental investigations of the process r — > \i ViVi 
would be informative in the determination of the physics beyond the SM and the free parameters 
existing in this model. 
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(a) 




Figure 1: One loop diagrams contribute to r — > fi Vi z/j, % = e,fi,r decay due to the neutral 
Higgs bosons h and A in the model III version of 2HDM. Solid lines represent leptons and 
neutrinos, curly (dashed) lines represent the virtual Z boson (h and A fields). 
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Figure 2: £§ TT dependence of the BR for real couplings and m h o = 85 GeV, m A o = 90 GeV. 
Here solid (dashed, small dashed, dotted, dash-dotted) line represents the case for €§ T/Ji = 
5 GeV (10, 15, 20, 25 GeV). 




Figure 3: The possible values of £ 



x T and the ratio je 



for m h o = 85 GeV, m A o = 90 GeV 



and the fixed values of the BR, BR = 10 4 (solid line) and BR = 10 6 (dashed line) 
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Figure 4: m h o dependence of the BR for the fixed values of Cn,tu 
and m^o = 90 GeV 
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Figure 5: The sin 6 TT dependence of the BR for m/jo = 85 GeV, m^o 
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